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Stellar evolution with rotation XII: Pre-supernova models 
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Abstract. We describe the latest developments of the Geneva stellar evolution code in order to model the pre- 
supernova evolution of rotating massive stars. Rotating and non-rotating stellar models at solar metallicity with 
masses equal to 12, 15, 20, 25, 40 and 60 Mq were computed from the ZAMS until the end of the core silicon 
burning phase. We took into account meridional circulation, secular shear instabilities, horizontal turbulence 
and dynamical shear instabilities. We find that dynamical shear instabilities mainly smoothen the sharp angular 
velocity gradients but do not transport angular momentum or chemical species over long distances. 
Most of the differences between the pre-supernova structures obtained from rotating and non-rotating stellar 
models have their origin in the effects of rotation during the core hydrogen and helium burning phases. The 
advanced stellar evolutionary stages appear too short in time to allow the rotational instabilities considered in 
this work to have a significant impact during the late stages. In particular the internal angular momentum does 
not change significantly during the advanced stages of the evolution. We can therefore have a good estimate of 
the final angular momentum at the end of the core helium burning phase. 

The effects of rotation on pre-supernova models are significant between 15 and 30 Mq . Indeed, rotation increases 
the core sizes (and the yields) by a factor ~ 1.5. Above 20 Mq, rotation may change the radius or colour of the 
supernova progenitors (blue instead of red supergiant) and the supernova type (Ibc instead of II). Rotation affects 
the lower mass limits for radiative core carbon burning, for iron core collapse and for black hole formation. For 
Wolf-Rayet stars (Af > 30 Mq), the pre-supernova structures are mostly affected by the intensities of the stellar 
winds and less by rotational mixing. 
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^ 1. Introduction stages. In Sect. 3 we present the stellar evolution in the 
r ""j Hertzsprung-Russell diagram and the lifetimes of the dif- 
Pg : Over the last years, the development of the Geneva evo- ferent burning stages. In Sect. 4 and 5 we discuss the evolu- 
^ lutionary code has allowed the study of rotating star evo- tion of rotation and internal structure respectively. Section 
■ - - lution from the ZAMS until the end of the core car- g describes the structure of the pre-supernova models, 
bon burning phase. Various checks of the validity of the Finally, in Sect. 7, we compare our results with the liter- 
rotating stellar models have been made. In particular, ature 
it has been shown that rotating models well reproduce 
the observed surface enrichments (|Heger fc Langenl2000t 
[Mcvnct & Macdcr 2000), the ratios of blue to red super- 
giant s in the Small Magellanic Cloud iJMaeder fc Mevnetl 
l200l|) . and the variations of the Wolf-Rayet (WR here- 
inaf ter) star populations as a function of the metallic- The comp uter model used here is th e same as the one de- 
ity l)Mevnet k MaedeJl2003(l. For all these features non- scribed in iMevnet fc MaedeJ l)2003|) except for the wind 
rotating models cannot reproduce observations. The goal anisotropy which here is not taken into account. The 
of this paper is to follow the evolution of these mod- model therefore includes secular shear and meridional 
els, which well reproduce the above observed features, circulation. Convective stability is determined by the 
during the pre-supernova evolution. Section 2 describes Schwarzschild criterion. Overshooting is only considered 
the modifications done in order to model the advanced for H- and He-burning cores with an overshooting param- 
eter, aovcr, of 0.1 Hp. The modifications made in order to 



2. Model physical ingredients 
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Raphael . Hirschi@obs . unige . ch below. 



R. Hirschi et al: Rotation: pre-supernova evolution 



2.1. Internal structure equations 

The inte rnal structure equations used are described in 
iMevnet fc Maeder (1997 ). These equations have been 
discretised acc ording to Sugimoto's prescription (see 
ISugimotdll97Cl|) in order to damp instabihties which de- 
velop during the advanced stages of stellar evolution. 
We note that the equations are still hydrostatic (no 
acceleration term) as in the pre-supernova models of 
[Limongi fc Chieffi (2003) . 



2.2. Nuclear reaction network 

The choice of the nuclear reaction network is a compromise 
between the number of chemical elements one wants to fol- 
low and the computational cost (CPU and memory). The 
network used for hydrogen (H) and h e lium ( He) burnings 
is the same as in iMevnet fc Maeden pOD.^ . For carbon 
(C), neon (Ne), oxygen (O) and silicon (Si) burnings, we 
chose to minimise the computational cost without losing 
accuracy for the energy production and the evolution of 
the abundance of the main elements. For this purpose, the 
chemical species followed during the advanced stages are 
a, 12C, 160, 20Ne, 24Mg, ^SSi, 32s, 36 Ar, 40Ca, ^4x1, 48Cr, 
^^Fe and ^^Ni. This network is usually called an a chain 
netwo rk. Note that |Timme s c t al. ( 20001) and Hix et alj 
l|1998(l show that even a network of seven elements is suf- 
ficient for this purpose. 

The system of equations describing the changes of the 
abundanc es by the nuclear react ions is resolved by the 
method of lArnett fc TrurarJ(|l969r) . This method has been 
chosen because it is very stable and rapid. It is therefore 
suitable to be included in an evolutionary code. We en- 
sured that we used small enough time steps to keep it very 
accurate. The use of small time steps ensures on top of it 
a good treatment of the interplay between nuclear burn- 
ing and diffusion since these two phenomena are treated 
separately (in a serial way) although they occur simulta- 
neously. 

The reactions r ates are taken from the NACRE 
iJAngulo et alJ Il999|) compilation for the experimen- 
tal reaction rates and from the NACRE website 
I http://pntpm.ulb.ac.be/nacre.htm \ for the theoretical 
ones. 

The nuclear energy production rates are derived from 
the individual reaction rates for C, Ne and 0-burning 
stages. During Si-burning, two quasi-e quilibrium groups 
form around ^^Si and ^^Ni respectivelv. iHix et alJ l)l998tl 
therefore only follow explicitly the reactions between ""^^Ti 
and "i^Cr and assume nuclear statistical equilibrium be- 
tween the other elements heavier than ^^gi. They choose 
the reaction between ^^Ti and ^*Cr because it is the bot- 
tleneck between the two quasi-equilibrium groups. We fol- 
lowed explicitly the abundance evolution of the 13 ele- 
ments cited above. Howeve r, for the energy production, 
we followed the method of iHix et al.l ()l998|) during Si- 
burning. We therefore only considered the reaction rate 



between 44x1 and "^^Cr and multiply them by the energy 
produced by the transformation of ^^Si into ^^Ni. 

2.3. Dynamical shear 

The criterion for stability against dynamical shear insta- 
bility is the Richardson criterion: 



Ri 
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N^ 1 
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where U is the horizontal velocity, z the vertical coordi- 
nate and N'^ the Brunt- Vaisala frequency: 
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where g is the gravity, S = —d\np/d\n.T)^,p, Hp is 
the pressure scale height, Vad = dlnT/dlnP)ad, V = 
dhiT/dlnP, V^ = dln/i/dlnP and ip = d\np/d\iifi)T,p. 

The critical value, Ric = 1/4, corresponds to the 
situation where the excess kinetic energy contained in 
the differentially rotating layers is equal to the work 
done against the restoring force of the density gradient 
(also called buoyancy force). It is therefore used by most 
authors as the limit for the occurr ence of t h e dyn amical 
shear. However, recent studies by ICanutol l)2002|) show 
that turbulence may occur as long as Ri < Ric ^ 1- This 
cri tical value is consistent w i th num erical simulations done 
by iBriJggen fc Hillebrandtl l|200lh where they find shear 
mixing for values of Ri greater than 1/4 (up to about 1.5). 

Different dynamical shear diffusion coefficients, D, can 
be found in the literature. Heger et al. (2000) use: 

Ri 
D = [imn{dinst,Hp]{l - max{— -, 0})]^/Tdy„ (3) 

nic 



where Tdyn = y^r-^/{Gmr) is the dynamical timescale and 
dinst the spatial extent of the unstable region, which is 
limited to one Hp. 



iBriJggen fc Hillebrandtl (|200l|) use another formula 
and they do numerical simulations to study the depen- 
dence of D on Ri. They find the following result: 



D = 



0.6 10^ 



Ri 



(4) 



2.3.1. The recipe 



The following dynamical shear coefficient is used, as sug- 
gested by J. -P. Zahn (priv. comm.): 



D = l.l 
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where r is the mean radius of the zone where the insta- 
bility occurs, AQ is the variation of Q over this zone and 
Ar is the extent of the zone. The zone is the reunion of 
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consecutive shells where Ri < Ric- This is valid if Pe > Ij 
where P^, the Peclet Number, is the ratio of cooling to 
dynamical timescale of a turbulent eddy. We calculated 
three fjni = 300kms~-^ 15 Mq models to see the impact 
of d ynamical shear and the importance of the value of 
Ric iJHirschi et al.ll2003a|) : one without dynamical shear, 
one with Ric = 1/4 and the last one with Ric = 1. See 
Sect. l4.1l for a discussion of the results. In the present grid 
of pre-supernova models, the dynamical shear is included 
with Ric = 1/4. 

2.3.2. Solberg-H0iland instability 

Solberg-H0iland stability criterion corresponds to the in- 
clusion of the effect of rotation (variation of centrifugal 
force) in the convective stability criterion. It is a combi- 
nation of the Ledo ux (or possibly Schwarzschild) and the 
Rayle igh criteria l|Maeder fc Mevneq l200(J: iHeeer et alJ 
1200011 . Both the dynamical shear and Solberg-H0iland in- 
stabilities occur in the case of a very large angular velocity 
decrease outwards (usual situation in stars, see Fig. EJ. 
Note that if there is a large increase outwards, dynami- 
cal shear instability occurs but not the Solberg-ll0iland 
instability. 

Both instabilities, shear instability and Solberg- 
Il0iland stability, occur on the dynamical timescale. We 
therefore expect them to have similar effects. The question 
is which instability sets in first? By comparing the stability 
criteria of the dynamical shear and of the Solberg-ll0iland 
instability: 

1/4 (drJ/dr)^r^ < N^ dynamical shear 

-2 n[2 n + (dSl/dr) r] < N^ Solberg - H0iland, 

where il is the angular velocity, r the radius and iV^ 
the Brunt- Vaisala frequency, it can be demonstrated that 
whenever a zone is unstable towards the Solbcrg-ll0iland 
instability, it is also unstable towards the dynamical shear 
instability. Indeed: 



because 



1/4 {dn/dT^T^ > -2 n[2 n + (df^/dr) r] 

1/4 (dr2/dr)2r2 + 2 r2[2 f7 + {dfl/di) r] = 
l/4[(dr2/dr)r-t-4r2]2 >0 

This means that the treatment of the dynamical shear 
instability alone is sufficient (since the timescales are sim- 
ilar). We therefore did not include explicitly the Solberg- 
H0iland instability in our model. 

2.4. Convection 

Convective diffusion replaces instantaneous convection 
from oxygen burning onwards because the mixing 
timescale becomes longer than the evolution timescale 
at that point. The numerical method used for this pur- 
pose is the method used for rotational diffusive mixing 



iJMevnet et alJl2004|) . The mixing length theory is used to 
derive the corresponding diffusion coefficient. Note that 
mult i-dimensional studies h ave been started on this sub- 
ject l(Bazan fc Arnettlll998l) . 

3. Hertzsprung-Russell (HR) diagram and 
lifetimes 

Stellar models of 12, 15, 20, 25, 40 and 60 Mq at so- 
lar metallicity, with initial rotational velocities of and 
300 kms~^ respectively have been computed. The value 
of the initial velocity corresponds to an average velocity 
of about 220kms~^ on the Main Sequence (MS) which 
is ver y close to the observed average value (see for in- 
stance Ekud3Eii3)- The calculations start at the ZAMS 
for the 12, 15, 20 and 25 Mq models and at the end 
of central He-burning for the 40 and 60 Mq models 
(for these models, we ta ke over the calculations done by 
iMevnet fc Maedeill2003r . The calculations reach the end 
of central Si-burning with models of rotating stars and the 
end of shell Si-burning with models of non-rotating stars. 
For the non-rotating 12 Mq star, Ne-burning starts at 
a fraction of a solar mass from the centre but does not 
reach the centre and the calculations stop there. For the 
rotating 12 Mq star, the model stops after 0-burning. 

The major characteristics of the models are sum- 
marised in Tables Q] and [3 In order to calculate lifetimes 
of the central burning stages, we take the start of a burn- 
ing stage when 0.003 in mass fraction of the main burning 
fuel is burnt. We consider that a burning stage is finished 
when the main fuel mass fraction drops below 10~^. The 
results would be the same if we had chosen 10^^ or 10"^. 
Neon burning is an exception because neon abundance 
does not drop significantly before the end of oxygen burn- 
ing. We therefore consider the end of Ne-burning when 
its abundance drops below 10^'^. Therfore the lifetimes 
for Ne-burning are to be considered as estimates. Other 
authors use the duration of the convective core as the life- 
time. We note that using the duration of convective cores 
as central burning lifetimes instead of threshold values of 
the central abundance of the main fuel would yield results 
very similar to those given in Table ^ for H, He, O and 
Si-burning stages. The core sizes are given at the end of 
central silicon-burning and at the last model calculated 
(which corresponds to a different evolutionary stage in 
the non-rotating and the rotating models as seen above). 
The inner limit of each core is the star centre. The outer 
limit is the point in mass where the sum of the mass frac- 
tion of the main burning products (helium for M^, carbon 
and oxygen for Mco, ^^Si-^-^Ti for Mgi and "^Cr-^eNi for 
Mpo) becomes less than 0.75 (superscript 75) or 0.50 (su- 
perscript 50). Another possibility to define the outer limit 
of a core is to consider the lagrangian mass where the 
mass fraction of the main fuel (helium for the CO cores) 
drops below 10^^. The CO cores thus obtained are given in 
Tables 1 and 2 (superscript 01). These limits are suitable 
for most masses (see Fig. I15II . However, for very massive 
stars (see Fig. I13|l . shell He-burning transforms most he- 
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Table 1. Initial properties and lifetimes of central burning stages of solar metallicity models. Also given are the 
total mass and the different core masses at the end of central silicon burning as well as at the last time step of our 
calculations. These last models correspond approximatively to the end of the first shell silicon burning for non-rotating 
models and slightly later then central silicon burning for rotating ones. All masses are in solar mass units. Lifetimes 
are in years with exponent in brackets (2.14 (-2)=2.14 10^^). Velocities are in kms~^. 



Initial model properties 


M'zAMS 


15 


15 


20 


20 


25 


25 


40 


40 


60 


60 


VZAMS 





300 





300 





300 





300 





300 


Lifetime of burning stages 


TH 


1.13 (7) 


1.43 (7) 


7.95 (6) 


1.01 (7) 


6.55 (6) 


7.97 (6) 


4.56 (6) 


5.53 (6) 


3.62 (6) 


4.30 (6) 


-THe 


1.34 (6) 


1.13 (6) 


8.75 (5) 


7.98 (5) 


6.85 (5) 


6.20 (5) 


4.83 (5) 


4.24 (5) 


3.85 (5) 


3.71 (5) 


TC 


3.92 (3) 


1.56 (3) 


9.56 (2) 


2.82 (2) 


3.17 (2) 


1.73 (2) 


4.17 (1) 


8.53 (1) 


5.19 (1) 


5.32 (1) 


TNe 


3.08 


0.359 


0.193 


8.81 (-2) 


0.882 


0.441 


4.45 (-2) 


6.74 (-2) 


4.04 (-2) 


4.15 (-2) 


TO 


2.43 


0.957 


0.476 


0.132 


0.318 


0.244 


5.98 (-2) 


0.176 


5.71 (-2) 


7.74 (-2) 


TSi 


2.14 (-2) 


8.74 (-3) 


9.52 (-3) 


2.73 (-3) 


3.34 (-3) 


2.15 (-3) 


1.93 (-3) 


2.08 (-3) 


1.95 (-3) 


2.42 (-3) 










End of central silicon 


burning 










Mtotsl 


13.232 


10.316 


15.694 


8.763 


16.002 


10.042 


13.967 


12.646 


14.524 


14.574 


Ml^ 


4.211 


5.677 


6.265 


8.654 


8.498 


10.042 


13.967 


12.646 


14.524 


14.574 


Mif^ 


2.441 


3.756 


4.134 


6.590 


6.272 


8.630 


12.699 


11.989 


13.891 


13.955 


M°c'o 


2.302 


3.325 


3.840 


5.864 


5.834 


7.339 


10.763 


9.453 


11.411 


11.506 


Mi? 


1.561 


2.036 


1.622 


2.245 


1.986 


2.345 


2.594 


2.212 


2.580 


2.448 


M|« 


1.105 


1.290 


1.110 


1.266 


1.271 


1.407 


1.464 


1.284 


1.458 


1.409 


Last model 


Ml," 


1.842 


2.050 


2.002 


2.244 


2.577 


2.894 


2.595 


2.868 


2.580 


2.448 


ME 


1.514 


1.300 


1.752 


1.260 


1.985 


1.405 


2.586 


1.286 


2.440 


1.409 



hum into carbon and oxygen and one could also consider 
that Mco includes the whole star. In that case we suggest 
another definition of Mco, which we name Mqq, defined 



by Mi"<t = MOIq + Jmoi 



Xqo dm, where Xco is the sum 

of ^^C and ^^O mass fractions. This definition gives an in- 
termediate value between M^q and the total actual mass 
of the star. 



3.1. Hertzsprung-Russell (HR) diagram 

The models calculated in the present work follow th e 
same tracks as the models from lMevnet fc Maeden ll2003|) . 
This is expected since the only difference between the two 
sets of models is the inclusion of dynamical shear in the 
present models. Here we concentrate on the 20 AIq mod- 
els. For that purpose, we also calculated 20 Mq mod- 
els with initia l rotati on velocities of 100 and 200 kms^^ 
l|2Q03b'). 



l|Hirschi et alJ l2003bl. The tracks of the 20 Mq models 
are presented in Figs. ^ 12] £^nd |21 Figure ^ shows the 
evolutionary tracks of the four different 20 M0 stars in 
the HR-diagram. The non-rotating model ends up as a 
red supergiant (RSG) like the model of oth er groups (see 
iHeerer fc LangeiibOOOl: ILimongi et al.ll2000|) . However, the 
rotating models show very interesting features. Although 
the 100 km s^^ model remains a RSG, the 200 km s^^ 
model undergoes a blue loop to finish as a yellow-red su- 
pergiant whereas the 300 km s^^ model ends up as a blue 
supergiant (BSG). Thus rotation may have a strong im- 



Table 2. Same as Table [H for the 12 Mq models. The 
non-rotating model starts Ne-burning off-centre and the 
burning never reaches the centre. The unburnt Ne-0 core, 
is given by Mnc-o- 



Initial model properties 


MZAMS 12 


12 


VZAMS 


300 


Lifetime of burning stages 


TH 1.56 (7) 


2.01 (7) 


-me 2.08 (6) 


1.58 (6) 


TC 6.47 (3) 


6.09 (3) 


TNe 


1.138 


TO 


4.346 


End of calculation 


Mtotai 11.524 


10.199 


MP 3.141 


3.877 


M'(^^ 1.803 


2.258 


Mg^o 1-723 


2.077 


Ml° 0.805 


1.340 


Mnc-o 0.096 


- 



pact on the nature of the supernova progenitor (red, blue 
supergiant or even Wolf-Rayet star) and thus on some 
observed characteristics of the supernova explosion. For 
instance the shock wave travel time through the envelope 
is proportional to the radius of the star. Since RSG radii 
are about hundred times BSG ones, this travel time may 
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Fig. 1. HR-diagram for 20 Mq models: solid, dashed, 

dotted-dashed and dotted lines correspond respectively to ^ig- 3- ^cff vs central helium mass fraction for 20 M© 

t;ini= 0, 100, 200 and 300 kms'^ We also indicate the models: sohd, dashed, dotted-dashed and dotted lines cor- 



position of the progenitor of SN1993J. 

Central Helium abundance 

^■^^ ! & HR -diagram for 20 M models 



300 km/s 



L(Lo) 



5.00 



4.61 '- 
4.67 




4.T2 

Teff(K} 



3:^r 



X(He)c 



Fig. 2. 3D HR diagram with central helium mass fraction 
as the third dimension for non-rotating and rotating 20 
Mq models. 



differ by two orders of magnitude depending on the initial 
rotational velocity. 

When does the star evolve back to the blue after a RSG 
phase and why ? Figure [3 is a 3D plot of the HR-diagram 
(in the plane) and the extra dimension represents the cen- 
tral helium mass fraction, Xc(^He). The extra dimension 
allows us to follow H, He and post He-burnings within a 
same diagram. Indeed, Xc(^He) increases during the main 
sequence, then decreases during He-burning and finally is 



respond respectively to Wini= 0, 100, 200 and 300 kms 



equal to zero during the post He-burning evolution. We 
can see that: 

— For the non-rotating model. He-burning starts when 
the star crosses the HR-diagram (LogTos ~ 4) and 
the star only reaches the RSG stage halfway through 
He-burning. Finally, the star luminosity rises during 
shell He-burning. 

— For the Ujni — 300 km s^^ model, the star is more lu- 
minous and becomes a RSG before He-burning igni- 
tion. These two factors favour higher mass loss rates 
and the star loses most of its hydrogen envelope before 
He-burning is finished. Thus the star evolves towards 
the zone of the HR diagram where homogeneous he- 
lium stars are found, i.e. in the blue part of the HR 
diagram. We can see that the star track still evolves 
during shell He-burning. 

Figure O is a projection of Fig. |2in the Log Toff ver- 
sus Xc('*He) plane. Although less intuitive than the 3D 
plot, it is more quantitative and still allows us to fol- 
low the various burning stages described above. Figure 
shows that all the rotating models become RSG before 
the beginning of the He-burning phase. The 100 km s^^ 
model luminosity is lower than for the 300 km s^^ model 
and therefore less mass is lost during He-burning and the 
burning ends before the hydrogen envelope is removed. 
The star therefore remains a RSG. The 200kms~^ model 
evolution is similar to the 300 km s""'^ model but the ex- 
tent of its blue loop is smaller. At the end of He-burning 
for the 200 km s^^ model, LogTcs — 4.28 and the star 
becomes redder before C-burning starts. 
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Initial Mass [M^] 
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Fig. 4. Burning lifetimes as a function of the initial mass 
and velocity. Solid and dotted lines correspond respec- 
tively to rotating and non-rotating models. Long-dashed 
and dotted-dashcd lines are used for rotating and non- 
rotating Ne-burning lifetimes to point out that they are 
to be considered as estimates (see text). 



Although the models discussed here are for solar metal- 
licity, one can note that the behaviours of the models with 
iJini between 200 and 300kms~^ are reminiscent of the 
evolution of the progenitor of SN1987A. Let us recall that 
this supernova had a blue progenitor which evolve d from 
a RSG stage (see e.g. the review bv lArnett et al.l ll989). 
In Fig. n we also indicate the position of the progenitor 
of S N 1993J. SN 1993J pr obably belongs to a binary sys- 
tem l|Podsiaklowski et al.M l993). Nevertheless it has com- 
mon points with our z;i„i =2GGkms~^ 20 Mq model: the 
star model and the progenitor of SN 1993J have approxi- 
mately the same metallicity, they have a similar position 
in the HR-diagram taking into account the uncertainties 
and they both have a small hydrogen rich envelope, mak- 
ing possible a change from type II to type lb some time 
after the explosion. 

3.2. Lifetimes 

The lifetimes are presented in Tables ^ and |21 and plotted 
in Fig. 21 We focus here our discussion on the effects of 
rotation on the lifetimes of the advanced burning phases. 
A discussion of the earlier ste llar evolutionary phases 
can be found in p revious papers IjMevnet fc Maedeill2003t 
iHeger et al.ll2000l) . For C-burning onwards, we have two 
patterns: 




M < 30Mq: Since the He-burning temperature is higher 
in rotating stars, the Sr ratio is smaller at the end of 



Log r [cm] 

Fig. 5. il variations as a function of the radius inside 15 
Mq models: the dashed line is a profile from a model with- 
out dynamical shear and the solid line from a model with 
dynamical shear and Ric — 1 during core 0-burning. The 
long and short arrows indicate the zones where i?j < 1 in 
the model without and with dynamical shear respectively. 
Note that the profiles do not differ significantly. 

He-burning and therefore the C-burning lifetimes are 
shorter. If C-burning is less important, less neon is 
produced and neon burning is also shorter. The trends 
for O- and Si-burnings are similar. 
M > 30Mq : The rotating stars become more rapidly WR 
stars and are more eroded by winds. The central tem- 
peratures for rotating models are therefore equal or 
even smaller than for non-rotating models. This leads 
to higher ^ ratios, longer C- and Ne-burnings phases. 

These two groups correspond to mass ranges where rota- 
tional mixing {M < 30Mq) or mass loss (M > 30Mq) 
dominates the other process. 



4. Rotation evolution 

4.1. Dynamical shear 

As said in Sect. 12.3.11 we calculated three Wini = 
300kms~^ 15 Mq models to see the impact of dynam- 
ical shear and the importance of the value of Ric'. one 
model without dynamical shear, one with Ric — 1/4 and 
the last one with Ric = 1. In Fig. |S1 the variation of the 
angular velocity, fi, as a function of the radius is shown 
inside 15 Mq stellar models in the core 0-burning phase. 
Arrows indicate the zones which are unstable against 
dynamical shear instability. These zones remain unsta- 
ble during the whole post core He-burning phase. Our 
simulations show that the characteristic timescale of the 



R. Hirschi et al: Rotation: pre-supernova evolution 



t/3 

"o 
o 

> 

Of) 

c 

(0 



a 

o 



M = 25 Mq, Vj, 



= 300 kms-', Z = 0.02, a„„^^=0.1 




last model 

Si end 

Si ignition 



ignition 
C ignition 

He ignition 
H end 

ZAMS 



_L_ 



10 

ni„ 



[Mo] 



r^ M = 25 Mg, Vi, 



= 300 kms-i. Z = 0,02, a„„.=0.1 

' ' over 




Fig. 6. Angular velocity as a function of the lagrangian Fig. 7. Local specific angular momentum profiles for tiie 
mass coordinate, rrir inside the 25 Mq model (uini= 300 25 Mq model (wini= 300 kms^^) at different evolutionary 
kms^^) at various evolutionary stages. stages. 



dynamical shear (oc R^/D) is always very short when 
using Eq. (O for the dynamical shear diffusion coefh- 
cient. Indeed, we obtain diffusion coefficients between f 0^^ 
and 10^^ cm^ sec~^. This is in general one or two or- 
de rs of magnitude larger than usi ng t he expressions give n 
bv iBruggen fc Hillebrand^ l)200l[l or iHeeer et all l|200(]() . 
However, the extent of the unstable zones is very small, 
a few thousandths of Mq. Therefore the shear mainly 
smoothens the sharp fi-gradients as can be seen in Fig. 
ISlbut does not transport angular momentum or chemical 
species over long distances. The general structure and the 
convective zones are similar between the model without 
dynamical shear and the one with dynamical shear. 

Concerning the Richardson criterion, there is no sig- 
nificant difference between the models using Ric = 1/4 
and Ric = 1. Except for the 15 Mq model discussed in 
this subsection, all the other models were computed with 
Ric = 1/4. 

4.2. Angular velocity, $7, and momentum evolution 

Figure ini shows the evolution of fl inside the 25 Mq model 
from the ZAMS until the end of the core Si-burning phase. 
The evolution of Q results from many different processes: 
convection enforces solid body rotation, contraction and 
expansion respectively increases and decreases fl in order 
to conserve angular momentum, shear (dynamical and sec- 
ular) erodes f2-gradients while meridional circulation may 
erode or build them up and finally mass loss may remove 
angular momentum from the surface. If during the core 
H-burning phase, all these processes may be important, 



from the end of the MS phase onwards, the evolution of ^ 
is mainly determined by convection, the local conservation 
of the angular momentum and, for the most massive stars 
only during the core He-burning phase, by mass loss. 

During the MS phase, H. decreases in the whole star. 
When the star becomes a red supergiant (RSG), fl at the 
surface decreases significantly due to the expansion of the 
outer layers. Note that the envelope is gradually lost by 
winds in the 25 Mq model. In the centre, O significantly 
increases when the core contracts and then the il profile 
fiattens due to convection, fi reaches values of the order of 
1 s~^ at the end of Si-burning. It never reaches the local 
break-up angular velocity limit, 0^, although, when local 
conservation holds, Qr/^c oc r~^/^. 

Figure[71shows the evolution of the specific angular mo- 
mentum, jj. — 2/3 n^r^, in the central region of a 25 Mq 
stellar model. The specific angular momentum remains 
constant under the effect of pure contraction or expansion, 
but varies when transport mechanisms are active. One sees 
that the transport processes remove angular momentum 
from the central regions. Most of the removal occurs dur- 
ing the core H-burning phase. Still some decrease occurs 
during the core He-burning phase, then the evolution is 
mostly governed by convection, which transports the an- 
gular momentum from the inner part of a convective zone 
to the outer part of the same convective zone. This pro- 
duces the teeth seen in Fig. The angular momentum of 
the star at the end of Si-burning is essentially the same 
as at the end of He-burning (by end of He-burning, we 
mean the time when the central helium mass fraction be- 
comes less th an lO"^). Th is result is very similar to the 
conclusions of iHeger et al.l |2000) on this issue. They find 
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that the angular momentum profile does not vary substan- 
tially after C-burning ignition (see Sect. 17.51 for a compar- 
ison). It means that we can estimate the pre-supernova 
angular momentum by looking at its value at the end of 
He-burning. We calculated, for the 25 Mq model, the 
angular momentum of its remnant (fixing the remnant 
mass to 3 Mq). We obtained Acm = 2.15 lO^Ogcm^ s"^ 
at the end of He-burning and £rcm = 1.63 10^°gcni^ s^^ 
at the end of Si-burning. This corresponds to a loss of 
only 24%. In comparison, the angular momentum is de- 
creased by a factor ~5 between the ZAMS and the end of 
He-burning. This shows the importance of correctly treat- 
ing the transport of angular momentum during the Main 
Sequence phase. 

5. Internal structure evolution 

5.1. Central evolution 

Figure |H1 (left) shows the tracks of the 15 and 60 Mq 
models throughout their evolution in the central tempera- 
ture versus central density plane (Log Tc-Log pc diagram) . 
Figure IHI (right) zooms in the advanced stages of the 12, 
20 and 40 AIq models. It is also very instructive to look at 
Kippcnhahn diagrams fFigs.[TTIand ll2(l in order to follow 
the evolution of the structure. Figure El helps understand 
the cause of the movements in the LogTc-Logpc diagram. 
We clearly identify two categories of stellar models: those 
whose evolution is mainly affected by mass loss (with an 
inferior mass limit of about 30 Mq), and those whose 
evolution is mainly affected by rotational mixing (see also 
Sect. 1213). We can see that for the 12, 15, 20 Mq models, 
the rotating tracks have a higher temperature and lower 
density due to bigger cores. The bigger cores arc due to 
the effect of mixing, which largely dominates the struc- 
tural effects of the centrifugal force. On the other hand, 
for the 40 and 60 Mq models, mass loss dominates mixing 
effects and the rotating model tracks in the Log Tc-Log pc 
plane are at the same level or below the non-rotating ones. 
In order to understand the evolutionary tracks in the 
LogTc-Logpc plane, we need to look at the different 
sources of energy at play. These are the nuclear energy, the 
neutrino and photon energy losses and the gravitational 
energy (linked to contraction and expansion). The differ- 
ent energy production rates at the star center are plotted 
in Fig. O as a function of the time left until core collapse. 
Going from the left to the right of Fig. O the evolution 
starts with H-burning where eh dominates. In response, a 
small expansion occurs (Cg negative and very small move- 
ment to lower densities in the 15 Mq model during H- 
burning in Fig. IS)). At the end of H-burning, the star 
contracts non-adiabatically (T ~ p^^"^, every further con- 
traction is also non-adiabatic). The contraction increases 
the central temperature. This happens very quickly and is 
seen in the sharp peak of Eg between H- and He-burnings. 
When the temperature is high enough, He-burning starts, 
SHe dominates and contraction is stopped. Note that dur- 
ing the H- and He-burning phases, most of the energy is 



transferred by radiation on thermal timescale. After He- 
burning, neutrino losses (e^ < 0) overtake photon losses. 
This accelerates the evolution because neutrinos escape 
freely. During burning stages, the nuclear energy produc- 
tion stops the contraction if e„„ci ^ ~£u (see C-burning 
for the rotating model) or even provoke an expansion 
when Enuci > —£v (most spectacular during Si-burning). 
Central density decreases when the central regions expand 
(see Fig. IHJ. Once the iron core is formed, there is no 
more nuclear energy available while neutrino losses are 
still present and the core collapses. 

Figure 1101 shows the variation of the energy production 
as a function of the mass fraction inside a 20 Mq stellar 
model at a stage during the shell C-burning phase. At the 
different burning shells, expansion occurs due to positive 
nuclear energy production. In the outer part, contraction 
and expansion are controlled by the photon luminosity and 
therefore by the opacities. In the inner regions, the energy 
produced either by the nuclear reactions or by contrac- 
tion is evacuated by the neutrinos. In the non-rotating 
star, partial ionisation of helium I in the outermost lay- 
ers produces a peak in the opacity (k y). This induces 
an expansion of the star (e^ < 0). The situation is differ- 
ent for the rotating model because it has lost most of its 
envelope and temperatures are higher than the ionisation 
transition zone. 

Numerically, it is important to note that the largest value 
for energy production rates corresponds to the nuclear 
one. Its maximum value is therefore used in order to de- 
termine the evolutionary time steps in our code. 

5.1.1. The fate of the 12 Mq models 

By looking at the track of the 12 Mq models in Fig.|Sl we 
can see that rotation has a noticeable effect on the post 
C-burning phases. Indeed, the non-rotating model starts 
Ne-burning off-centre and the burning never reaches the 
centre. The unburnt Ne-0 core, Mno-o is equal to 0.096 
(see TableEJ- On the other hand, the rotating model starts 
Ne and 0-burnings in the centre. This can be seen in Fig. 
ITTI (top). The computation of the 12 Mq models were 
stopped during the Ne/0-burning phase. To explore their 
further evolution, on e can use the rnass limits for the Ne- 
cores, Mnb, given bv iNomota (|l9843: 

— AfiMc — 1-46 Mq is the lower limit for neon ignition in 
the centre. 

— Mnb = 1-42 Mq is the lower limit for off-centre 
neon ignition where the subsequent neon burning front 
reaches the centre. 

— Mnc = 1-37 is the lower limit for off-centre neon ig- 
nition. In the mass range between 1.37 and 1.42 Mq, 
neon burning never reaches the centre. 

Our rotating 12 Mq model has Mno ^ 1-6 well above 
the lower mass limit for neon ignition in the centre. In 
the non-rotating model the Ne-core mass is around 1.4 
Mq (more quantitatively carbon mass fraction decreases 
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Fig. 9. Log of the energy production rate per unit mass at the star center as a function of the time left until core 
collapse for the non-rotating (left) and rotating (right) 20 Mq models. Nuclear energy production rates during H- 
and He-burnings are shown in dotted (eh) and dashed (shc) lines respectively. The solid line corresponds to the 
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from 0.01 to 0.001 between 1.45 and 1.34 Mq). As ex- 
pected from the mass limits above, in this model, neon 
burning, which st arts off-centre, will probably reach the 
centre. Then fsee iNomoto fc Hashimotol ll988\ Sect. 3.2: 
fate of stars with 10 M(7) < Mms < 13 Mr^, and references 



^0 







therein), electron capture will help the star to collapse 
making the neon/oxygen burning explosive and possibly 
ejecting the H and He-rich layers. Note that in our mod- 
els we only follow multiple-a elements. We did not follow 
the evolution of the electron mole number, Y^., or of neu- 
tron excess, 77, neither include Coulomb corrections. Let 
us recall that the electron mole number, Yg ~ J2i^i^i^ 
and the neutron excess, 77 — J2ii-^i ~ ^iWi^ are linked 
by the following relation: Yg = (1 — rj)/2 {Ni, Zi and Yi 
are respectively the number of neutron s, protons and the 
number abundance of element i; Yi = Xi/Ai, where Xi 
and Ai are the mass fraction and the mass number of ele- 
ment i). Therefore the electron mole number, Yg, is always 
equal to 0.5""^. Lower values of Yg (due to electron captures) 
and the inclusion of Coulomb corrections in the equation 
of state have an impact in this context. Electron cap- 
tures remove electrons. This decreases the electron pres- 
sure and facilitates the collapse. Coulomb corrections gen- 
erally act to decrease the iron core mass by about 0.1 Mq 
{ Wooslev et al.ll2002l and references therein) . These omis- 
sions can be the cause of the failure of our models to follow 



^ The mass limits given bv lNomotd l|l98j) were also obtained 
from calculations with Y. = 0.5. 



the evolution of the 12 Mq models further. These two ef- 
fects however do not affect significantly the evolution of 
more massive stars before the shell Si-burning phase. 

5.2. Kippenhahn diagrams 

Fiefures ITTI II 21 show the Kippenhahn diagrams for the dif- 
ferent models. The y-axis represents the mass coordinate 
and the x-axis the time left until core collapse. The black 
zones represent convective zones. Since our calculations 
have not reached core collapse yet, we estimate that there 
is 10^^ yr between the last model and the collapse. This 
value has no significant influence since it is only a small 
additive constant. The graph is built by drawing vertical 
lines at each time step where the star is convective. this 
discrete construction shows its weakness at the right edge 
of each diagram and during shell He-burning where time 
steps are too distant from each other to cover the surface 
properly. The abbreviations of the various burning stages 
are written below the graph at the time corresponding to 
the central burning stages. 

We can see the effect of the blue loops 
llMevnet fc Maeded 12003^ in the 12 Mq models on 
the external convective zone during the core He-burning 
phase. The blueward motion reduces the external con- 
vective zone or even suppresses it. We also note the 
complex succession of the different convective zones 
between central O and Si-burnings (for instance in the 
non-rotating 15 Mq model). The difference between 
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non-rotating and rotating models is striking in the 20 
and 25 Mq models. We can see that small convective 
zones above the central H-burning core disappear in 
rotating models. Also visible is the loss of the hydrogen 
rich envelope in the rotating models. On the other hand 
non-rotating and rotating 40 and 60 Mq models all have 
very similar convective zones history after He-burning. 



5.2.1. Convection during core C-burning? 

Recent calculations (jHeger et al.l l2000|) show that non- 
rotating stars with masses less than about 22 Mq have a 
convective central C-burning core while heavier stars have 
a radiative one. Our non-rotating models agree with this. 
What about models of rotating stars? Figure [TTl fbottom) 
shows the Kippenhahn diagrams for the non-rotating and 
rotating 20 Mq models. We can see that the rotating 
model has a radiative core during central C-burning. It 
is due to the fact that the nuclear energy production rate 
ec does not overtake — e^ (see Fig. Alright) and therefore 
the central entropy does not increase enough to create a 
convective zone. This behaviour results from the bigger 
He-cores formed in rotating models. Bigger cores imply 
higher central temperatures during the core He-burning 
phase and higher central temperatures imply lower carbon 
content at the end of the He-burning phase. Thus less fuel 
is available for the core C-burning phase which does not 
succeed to develop a convective core. The same explana- 
tion works for more massive (rotating or non-rotating) 
stars. Thus the upper mass limit for a convective core 
during the C-burning phase is lowered by rotation, pass- 
ing from about 22 Mq to a value inferior to 20 Mq when 
the initial velocity increases from to 300 km s^^. 

5.3. Abundances evolution 

Figures ^^ and El show the evolution of the abundances 
inside the non-rotating (left) and rotating (right) 20 Mq 
models at the end of each central burning episode. A the 
end of H-burning, we notice the smoother profiles in the 
rotating model, consequence of the rotational mixing. At 
the end of He-burning, we can already see the difference 
in core sizes and total mass. We also notice the lower C/0 
ratio for rotating models. At the end of 0-burning, we can 
see that the rotating model produces much more oxygen 
compared to the non-rotating model (about a factor two). 
At the end of Si-burning, the iron and Si-cores are slightly 
bigger in the rotating model (see also Table ^. The yields 
of oxygen are therefore expected to increase significantly 
with rotation. This will be discussed in an forthcoming 
article. 



M = 60 Mq, Vi„i=300 kms-i, Z = 0.02, a^^^^--0.1 



CD 



CD 



X 





1 ' 


-''!■' 


. 1 . 1 


' ' ' J-, 


1 

1 


1 , 


: 


: ^--i 


160 






" / 


J 
! ! 

1 
1 i 


( 




V" 


•He 


; 


- 


1/ 





— 






r/, 


,■ .'1 






I 






^ 


V 


1 '\ i 
' ' 1 


i I — , 
li 
'.il 








L 


- 


™Si 


1 


;'/-\: 1 










1 ^ 




azS 








" 


^1 












; 


1 fj 




< 1 < 1 




1 


1 , 


- 



Fig. 13. Variations of the abundance (in mass fraction) 
as a function of the lagrangian mass coordinate, mr, at 
the end of central Si-burning for the rotating 60 Mq. 
Note that the ^''Ti abundance (dotted-long dashed line) 
is enhanced by a factor 1 000 for display purposes. 
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Fig. 16. Core masses as a function of the initial mass and 
velocity at the end of core Si-burning. 



6. Pre-supernova models 

6.1. Core masses 

Figure El shows the core masses (Tables H and |21l as a 
function of initial mass for non-rotating (dotted lines) 



and rotating (solid lines) models. Since rotation increases 
mass loss, the final mass, Mfi„ai, of rotating models is al- 
ways smaller than that of non-rotating ones. Note that 
for very massive stars (M > 60Mq) mass loss during the 
WR phase is proportional to the actual mass of the star. 



14 



R. Hirschi et al: Rotation: pre-supernova evolution 



M = 20 M„, i;,„, = 000 km s-', Z = 0.02, a„„=0.1 



■^He 



end H 



o 



(D 

X 



X 10-= 



O 



X 



X 10-' 



■*He 






/ 




1 m end H 


14N 

- 
- 
- 


J^--^^:^""" : 


/ 


- 160 ! 

: 1 


1 / : 


) 


: -c /■ 




< < < 1 < < < 1 < < < 1 < < < 1 < < < 1 < < < 1 < < 



6 B 10 



1-1 16 Ifl 20 



6 8 10 



1-1 16 IB 30 



M = 20 Mg, i;i„i = kms-', Z = 0.02, a^,„=0.1 



m, [M^: 




160 

ISC 



o 



end He 



^Ne 



X - 20Ne 



X 10-= 



■*Mg 




m, [MJ 
M = 20 M„, V. . = krns-i, Z = 0.02, a,„^ =0,1 



m, [M^: 




M = 20 M„, i;,„, = 300 kms-', Z = 0,02, a„„=0,l 



me 



o 



u , 



120 

end C 



iHe 



X 10-' r '-'1 



m, [Mj 



m, [Mj 
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This produces a convergence of the final masses (see for 
instance Meynet & Maeder 2004). We can also see a gen- 
eral difference between the effects of rotation below and 
above 30 Mq. For M < 30 Mq, rotation significantly in- 
creases the core masses due to mixing. For M > 30 Mq , 
rotation makes the star enter at an earlier stage into the 
WR phase. The rotating star spends therefore a longer 
time in this phase characterised by heavy mass loss rates. 
This results in smaller cores at the pre-supernova stage. 
We can see on Fig.Elthat the difference between rotating 
and non-rotating models is the largest between 15 and 25 
Mq. 

Concerning the initial mass dependence, one can make 
the following remarks: 

-^^finai: There is no simple relation between the final mass 
and the initial one. The important point is that a final 
mass between 10 and 15 Mq can correspond to any 
star with an initial mass between 15 and 60 Mq. 

Ma and Mqo ■ The core masses increase significantly with 
the initial mass. For very massive stars, these core 
masses are limited by the very important mass loss 
rates undergone by these stars: typically Ma is equal 
to the final mass for M > 20 Mq for rotating models 
and for M > 40Mq for the non-rotating ones. The 
mass of the carbon-oxygen core is also limited by the 
mass loss rates for Af > 40 Mq for both rotating and 
non-rotating models. 

Msi (at the end of central Si-burning): For rotating 

models, A/si oscillates between 2 and 2.5 Mq. For 
non-rotating models, the mass increases regularly 
between 15 Mq (Msi ~ 1.56Mq) and 40 Mq 
(Msi ~ 2.6Mq) and stays constant for higher masses 
(due to mass loss). 

Mpc (at the end of central Si-burning): Follows the same 
trend as Afsi. 

6.1.1. Final iron core masses 

For non-rotating models, the masses of the iron core Afpc 
in the last computed model (end of shell Si-burning) 
are very close (within 8%) to the silicon core masses, 
Msi, at the end of central Si-burning. This occurs be- 
cause the extent of shell Si-burning is limited by the en- 
tropy increase produced by the second episode of shell 
0-burning. Therefore even though our rotating models 
have not reached core collapse, we can have an estimate 
of the final iron core mass by taking the value of M^i at 
the end of central Si-burning. In this way, we obtain iron 
core masses for rotating models between 2 and 2.5 Mq. 
For non-rotating models, the mass is between 1.56 Mq 
and 2.6 Mq. Rotating models have therefore more mas- 
sive iron cores and we expect the lower mass limit for black 
hole (BH) formation to decrease with rotation. 

As said in Sect. 15.111 about the fate of the 12 Mq, we 
did not follow the evolution of the electron mole num- 
ber. Ye, neither include Coulomb corrections. Coulomb 
corrections generally act to decrease the iron core mass 



by about 0.1 Mq JWooslev et al.l l2002l and references 
therein). Electron captures during Si-burning increases 
neutron excess and also reduces the electron pressure and 
this ( with photodisint egration) will allow the core to col- 
lapse (jWooslev et all 12002'). It is therefore possible that 
some of our models should collapse before shell Si-burning 
occurs. Taking this argument into account and the fact 
that we used Schwarzschild criterion for convection, we 
have to consider the value of Msi at the end of core Si- 
burning as an upper limit for the final iron core mass. 

6.2. Internal structure 

As well as the chemical composition (abundance profiles 
and core masses) of the pre-supernova star, other param- 
eters, like the density profile, the neutron excess (not fol- 
lowed in our calculations), the entropy and the total radius 
of the star, play an important role in the supernova explo- 
sion. Figure lTTI shows the density, temperature, radius and 
pressure variations as a function of the lagrangian mass co- 
ordinate at the end of the core Si-burning phase. Since the 
rotating star has lost its envelope, this truly affects the pa- 
rameters towards the surface of the star. The radius of the 
star (BSG) is about one percent that of the non-rotating 
star (RSG). As said above this modifies strongly the su- 
pernova explosion. We also see that temperature, density 
and pressure profiles are flatter in the interior of rotating 
models due to the bigger core sizes. 

7. Comparison with the literature 

In this section, we compare our results (HMM here- 
inafte r) with four other recent papers : [ Lirnongi et alJ 
('200 01) (LSC he reinafter) , "moieiieraOjliooJiWHWr 
Rau scher et all (|2b02) (RHWl and iHeeer et all l|2000h 
(HLW). Before we start the comparison, we need to men- 
tion which physical ingredients (treatment of convection, 
^'^C{a,j)'^^0 reaction rate, . . . ) they use: 

— LSC use Schwarzschild criterion for convection with- 
out overshooting (except for core He-burning for which 
semiconvection and an induced overshooting are taken 
into account). For ^^C (a,7)''^^0, they use the rate of 
ICaughlan et all l)l985|) (CF85). Mass loss is not in- 
cluded. 

— WHW use Ledoux criterion for convection with semi- 
convection. They use a relatively large diffusion co- 
efficient for modeling semiconvection. Moreover non- 
convective zones immediately adjacent to convective 
regions are slowly mixed on the order of a radiation 
diffusion time scale to approximately allow for the ef- 
fects of conve ctive overshoot. For _|^C(a,7)^^0, they 
use the rate of lCaughlan fc Fowled l|l988l) (CF88) mul- 
tiplied by 1.7. 

— RHW use Ledoux criterion for convection with semi- 
convection. They use the same method as WHW for 
semiconvection. For ^^C(a,7)^^0, they use the rate of 
iBuchmanpi 11996) (BU96) multipHed by 1.2. 
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Fig. 17. Profiles of the radius, r, density, p, temperature, T and pressure P at the end of core Si-burning for the 
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— HLW use Ledoux criterion for convection with semi- 
convection using a small diffusion coefhcient (about 
one percent of WHW's coefficient) and without over- 
shooting. For ^ ^ C(q;,7 )^^0, they use a rate close to 
ICaughlan erall l)l985|) (CF85). They present models 
with and without rotation. 

— In this paper (HMM), we used Schwarzschild crite- 
rion for convection with overshooting for core H and 
He-burnings. For ^ '^Cia^^Y^O. we used the rate of 
lAngulo et ail (|l 99fl^ (NACRE). 

7.1. HR diagram 

We remark that the present evolutionary tracks (as well 
as those from LSC) do not decrease in luminosity when 
they cross the Hertzsprung gap. This is in contrast with 
the tracks from HLW which present a significant decrease 
in luminosity when they evolve from the MS phase to the 
RSG phase. Models computed with the present code but 
using the Ledoux criterion for convection (without semi- 
convection) present a very similar behaviour to those of 
HLW. Thus the difference between the two sets of models 
mainly results from the different criterion used for convec- 
tion. 



7.2. Lifetimes 

We can compare the lifetimes of the non-rotating 15, 
20, 25 Mq models with recent calculations presented in 
WHW and LSC. The comparison is shown in Table |21 As 
said earlier, LSC use Schwarzschild criterion with over- 
shooting only for He-burning. WHW use Ledoux criterion 



with a very efficient semiconvection and allow for some 
overshoot. Despite important differences in the treatment 
of convection, all the models give very similar H-burning 
lifetimes which differ by less than 10%. For the He- 
burning lifetimes, during which the convective core grows 
in mass, one can expect that the results will be signifi- 
cantly different depending on which convection criterion 
is used. This is indeed the case. Inspecting Table 3, one 
sees that our results are shorter by 30-50% with respect to 
those of WHW. In contrast when the Schwarzschild crite- 
rion is used with some overshooting as in LSC, the results 
are very similar (differences inferior to six percents). In the 
advanced stages one sees that the lifetimes obtained by the 
different groups are of the same order of magnitude. Let 
us note that the definition of the duration of the nuclear 
burning stages may differ between the various authors and 
this tends to enhance the scatter of the results. Keeping 
in mind this source of difference and the fact that the life- 
times vary by eight or nine orders of magnitude between 
the H-burning and the Si-burning phases, the agreement 
between the various authors appears remarkable. 



7.3. Kippeniiahn diagrams and convection during 
central C-burning 

Our Kippenhahn diagrams for th e non-rotating models 
are in good agreement with those o flRauscher et alJ l|2002r) 
except that in our model the carbon and oxygen shells do 
not merge for the 20 Mq. The only noticeable difference 
between the structures in the advanced phase obtained 
in the present work and those obtained by LSC is that 
their 20 Mq model does not have a central convective 
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Table 3. Lifetimes of central burning stages of solar metallicity models. Lifetimes are in years with exponent in 
brackets (2.14 (-2)^2.14 lO^^). 



MzAMS 


15 (HMM) 


15 (WHW) 


15 (LSC) 


20 (HMM) 


20 (WHW) 


20 (LSC) 


25 (HMM) 


25 (WHW) 


25 (LSC) 


TH 


1.13 (7) 


1.11 (7) 


1.07 (7) 


7.95 (6) 


8.13 (6) 


7.48 (6) 


6.55 (6) 


6.70 (6) 


5.93 (6) 


THe 


1.34 (6) 


1.97 (6) 


1.4 (6) 


8.75 (5) 


1.17 (6) 


9.3 (5) 


6.85 (5) 


8.39 (5) 


6.8 (5) 


rc 


3.92 (3) 


2.03 (3) 


2.6 (3) 


9.56 (2) 


9.76 (2) 


1.45 (3) 


3.17 (2) 


5.22 (2) 


9.7 (2) 


TNe 


3.08 


0.732 


2.00 


0.193 


0.599 


1.46 


0.882 


0.891 


0.77 


TO 


2.43 


2.58 


2.43 


0.476 


1.25 


0.72 


0.318 


0.402 


0.33 


TSi 


2.14 (-2) 


5.01 (-2) 


2.14 (-2) 


9.52 (-3) 


3.15 (-2) 


3.50 (-3) 


3.34 (-3) 


2.01 (-3) 


3.41 (-3) 



Table 4. Final core masses at the pre-supernova stage for 
different models of non-rotating stars at solar metallicity. 



Core masses vs initial mass 



M; 



ZAMS 15 (HMM) 15 (RHW) 15 (HLW) 15 (LSC) 



M:otal 


13.232 


12.612 


13.55 


15 


M«i 


4.168 


4.163 


3.82 


4.10 


Mg^o 


2.302 


2.819 


1.77 


2.39 


Ml? 


1.842 


1.808 


- 


- 


M^2 


1.514 


1.452 


1.33 


1.429 


MzAMS 


20 (HMM) 


20 (RHW) 


20 (HLW) 


20 (LSC) 


Mtotal 


15.694 


14.740 


16.31 


20 


Mf 


6.208 


6.131 


5.68 


5.94 


M°c'o 


3.840 


4.508 


2.31 


3.44 


Ml? 


2.002 


1.601 


- 


- 


Mf 


1.752 


1.461 


1.64 


1.552 


MzAMS 


25 (HMM) 


25 (RHW) 


25 (HLW) 


25 (LSC) 


Mtotal 


16.002 


13.079 


18.72 


25 


M°i 


8.434 


8.317 


7.86 


8.01 


M'c'o 


5.834 


6.498 


3.11 


4.90 


Ml? 


2.577 


2.121 


- 


- 


Mf 


1.985 


1.619 


1.36 


1.527 



core during C-burning. This can be expl ained by the fact 
that t hey use the ^^C{a,"f)^^0 rate from lCaughlan et al.1 
||198,'t|^ . This rate is larger than the NACRE rate we use in 
our models (see Fig. l21|) and therefore more ^^C is burnt 
during He-burning. 

7.4. Core masses 

7.4.1. Non-rotating models 

In Table 01 the final masses and the core masses at the 
pre-supernoya stage are given for different 15, 20 and 25 
Mq stellar models. The second column corresponds to the 
pr esent non-rotat i ng mo dels, the third show s the results 
of Raus cher et al.l Ij2002l) , the f ourth those o f Heeer et al. 
(2000) and the fifth those of iLimongi'eFal. (2000) . In 
Fig. ^1 we see that for Ma, the results are very similar 
(within 5%) between our models and those of LSC and 
RHW. This can be understood by the similar outcome of 
the convection treatment. HLW use a small diffusion co- 
efficient for semiconvection and logically obtain slightly 
smaller helium cores. 

The differences between the mass of the CO cores are 
much greater. Let us recall here that the size of this core 
depends a lot on the convective criterion and also on the 
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Fig. 18. Variation of the He core masses, Ma (light lines), 
and of the CO core masses, Mco (heavy lines), at the 
pre-supernova stage in different initial mass models. Only 
non-rotating models are shown. The different types of 
line correspond to results obtained by different groups: 
HMM labels our result s, W86 those of ^yooslev fc Weaven 
ill 9861). RHW those oflRanischer et aLl ll2002f) LSC those 
of iLimongj et al.l I 200(T() and HLW those of iHeger et alJ 
(l200(t. 



rate of the ^^C(a,7)^®0 reaction. This reaction becomes 
one of the main source of energy at the end of the core 
He-burning phase. A faster rate implies smaller central 
temperatures and thus increases the He-burning lifetime. 
This in turn will produce larger CO cores iJLanger 1991() 
with a smaller fraction of ^^C. Figure l?T1 shows the rates 
used by various authors divided by the NACRE rate for 
the temperature range of interest. 

Since HLW and LSC use the same rate for this reac- 
tion, most of the difference between the mass of the CO 
cores must have its origin in the different treatment of 
convection. One notes also that LSC still have slightly 
smaller cores than us even though they added some 
semiconvection and use the CF85 rate for ^^C(a,7)^^0 
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Fig. 19. Miron (and Afsi) as a function of Ma for non- 
rotating models from different authors (see Table 0J. The 
labels are the same as in Fig. ^| The light lines show the 
variation of A'/si, the heavy lines those of Mjion- 
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Fig. 21. Comparison of different ^^C(a,7)"'^^0 reaction 
rates accordi ng to various authors : CF85 labels the 
rate given bv ICaughlan et all ||198,'t() . CF88, the rate of 
ICaughlan fc Fowled l|l988l). BU9 6 the rate of'B uchmani] 
(!:1996) and K02 the one of lKunz"et al. ( 2002,) AH the rates 
are no rmalised to the rate given by NACRE IjAngulo et alJ 
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Fig. 20. Miron (and Msi) as a function of Mco for non- 
rotating models from different authors (see Table ^. The 
labels are the same as in Fig. El The light lines show the 
variation of Msi, the heavy lines those of A/iron- 



which is greater than the one adopted in our models. 
RHW, although they had slightly smaller Ma , have larger 
Mco- This can be explained in part by the use of the 



rate BU96xl.2 for ^^C(a,7)^^0 which is larger than the 
NACRE rate at the end of He-burning (see Fig.l?T|l. 

In Figsll9landl20lthe Si and iron core masses obtained 
at the end of shell Si-burning are plotted as a function 
of Ma and Mco respectively. The present results are well 
in the range of values obtained by different authors for 
Ma <6Mq (Mini < 20 Mq). Above t his mass range, our 
result s are in agreement with those of lWooslev fc Weaven 
l|198fTt) and significantly above the results obtained more 
recently by the other groups. As discussed in Sect. 16. I'll 
we did not follow the evolution of the neutron excess or 
include Coulomb corrections. This does not affect our re- 
sults until the end of core Si-burning but may affect the 
results plotted in Figs. 1191 and 1201 obtained at the end of 
the shell Si-burning. In this last case, the present results 
have to be considered as upper limits. This might be part 
of the explanation why our iron core masses appear to be 
systematically greater than tho se obtained in recent c alcu- 
lations. However one notes that lWooslev et al.l l 2002l) give 
a Chandrasekhar mass (lower mass limit for collapse) of 
1.79 Mq for the 25 Mq model which is large compared to 
the iron core we obtain at the end of core Si-burning, im- 
plying that our 25 Mq (Ma = 8.4 Mq) model may experi- 
ence shell Si-burning before collapsing. Thus if we cannot 
discard that the final iron core masses are overestimated 
due to the above reason, they may also be greater than 
the masses obtained by other groups for other reasons. In 
this context it is interesting to compare the masses of the 
Si-burning core. The Si-cores are created by 0-burning 
before Si-burning (except possibly a small fraction due to 
an additional shell 0-burning during Si-burning). Their 
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sizes are thus not dependent on the neutron excess or the 
Chandrasekhar mass. Looking at Fig. I19landl20l where our 
Si-core masses are compared to those obtained by RHW, 
we see that our core masses are systematically larger. In 
that case the difference cannot be attributed to the ne- 
glect in our models of the electron capture reactions and 
of the Coulomb corrections. Our bigger cores result from 
the different prescription we used for convection in our 
models. Thus it is possible that the bigger iron cores we 
have obtained are due, at least in part, to the prescriptions 
we used for convection. 



7.4.2. Rotating models 

We can also compare core masses of th e rotating 15, 20, 25 
Mq models with recent calculations bv lHeger et alJ l|200ni) 
(HLW). For A/fc, we use M^^ . As discussed in section 
16.1.11 this assumes that shell Si-burning occurs before the 
collapse and our value has to be considered as an upper 
limit. 

The comparison is shown in TablcEl "F..B" models are 
the models with the same initial rotational velocity and 
inclusion of the /z-gradients inhibiting effects on rotational 
mixing. These are the models which should give approxi- 
mately the same results as us if uncertainties concerning 
the treatment of convection and particular reaction rates 
were small. We also show the "E.." models with a lower 
initial rotational velocity but without the /i-gradients in- 
hibiting effects. 

One can see by comparing the results from the two 
models of HLW, the great dependence of the core masses 
on the treatment of the /x-gradient inhibiting effect. The 
more efffcient the rotational mixing (or less strong are the 
inhibiting effects of the /z-gradients) , the greater the core 
masses. Compared to the results obtained by HLW, one 
sees that our core masses are significantly greater. This 
essentially results from two facts: first the effects of rota- 
tion are included in our models in a different way than 
in the models by HLW. In particular in our models, the 
treatment of rotational mixing includes the inhibiting ef- 
fect of /x-gradients without any ad hoc parameters, and 
the transport of the angular momentum by the meridional 



M = 20 Mq, Z = 0.02, a„„^^=0.1 



circulation is properly accounted for by an advective term 
i Maeder fc Zahnll998l) . Secondly, the present stellar mod- 



els were computed with the Schwarzschild criterion with 
a moderate overshooting while the models of HLW were 
computed with the Ledoux criterion with semiconvection 
using a small diffusion coefficient and without overshoot- 
ing. 



7.5. Final angular momentum 

Long soft gamma-ray bursts (GRBs) were recently con- 
nected with SNe (see iMatheson 20031 for example). One 
scenario fo r GRB p r oduct ion is the coUapsar mechanism 
devised bv lWooslevI pOQ,*^ . In this mechanism, a star col- 
lapses into a black hole and an accretion disk due to the 
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200 km s-i 




m, [MJ 

Fig. 22. Comparison of the final local specific angular mo- 
mentum profiles for different 20 Mq models. Models with 
different initial velocities, Uini= 100, 200 and 300 km s^^ 
are drawn with dotted, solid and thick solid lines respec- 
tively. We can see the convergence of the final of the final 
angular momentum of the core above Uini= 200kms~^. 

high angular momentum of the core. Accretion from the 
disk onto the central black hole produces bi-polar jets. 
These jets can only reach the surface of the star (and 
be detected) if the star loses its hydrogen rich envelope 
before the collapse. WR stars are therefore good candi- 
dates for coUapsar progenitors since they lose their hy- 
drogen rich envelope during the pre-SN evolution. The 
question to answer is whether the core of WR stars con- 
tains enough angular momentum at the pre-SN stage (the 
specific angular momentum, j, of the material just out- 
side the core must be larger than 10^^ cm^s~^). So far 
only Heger and co-workers have obtained values for the 
angular momentum of the cores of m assive stars at the 
pre-SN stage (Heger et al. 200d l2003|) . The physical in- 
gredients of their model have been given in Sect. 17.4.21 
The comparison between our models and theirs shows that 
the size of the various cores depends significantly on the 
treatment of both convection and rotation. The evolution 
of angular velo city and angular momentum in the models 
of iHeger et al.i (2000. 2003) , is presented in Hcger et alJ 
l|200(i HLWOO hereinafter; in Figs. 8 and 9). The evolu- 
tion of angular ve locity and mo mentum in our models is 
described in Mevnet fc Maeded (2000) and in Sect.Ol 

HLWOO show with their Fig. 9 the convergence of the 
final angular momentum of the core for a wide range 
of initial angular momentum. The dependence of the fi- 
nal angular momentum on the initial one for our models 
is displayed in Fig. EH Models with Vini^ 100 and 200 
kms~^ have been computed until the end of 0-burning. 
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Table 5. Final core masses at the pre-supernova stage for different models of rotating stars at solar metallicity. Note 
that we use Msi at the end of central Si-burning for the value of Mpe as discussed in the text. 



MzAMS 


15 


F15B (HLW) 


E15 (HLW) 


20 


F20B (HLW) 


E20 (HLW) 


25 


E25 (HLW) 


Mtotal 


10.316 


12.89 


10.86 


8.763 


14.76 


11.00 


10.042 


5.45 


M«i 


5.604 


3.88 


5.10 


8.567 


5.99 


7.71 


10.042 


5.45 


Mg^o 


3.325 


2.01 


3.40 


5.864 


2.75 


5.01 


7.339 


4.07 


M|» 


2.036 


1.38 


1.46 


2.245 


1.36 


1.73 


2.345 


1.69 
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Fig. 23. Comparison of the final local specific angular 
momentum profiles for different 20 Mq models all with 
the same initial rotational velocity, Vini= 200 kms^^. The 
thick solid l ine (HMM) corresp onds to our model. The 
models from lHeeer et al.l ()2000|) are drawn with a dashed 
line for model E20 (no /i-barrier) and with a dotted- 
dashed line for model E20 B (jit-barrier with f „ — 0.05). 
Finally, model m20b5 from IHecfer et al.l 11200,1) including 
the effect of the magnetic fields according to lSpruit 1.2002) 
is drawn with the dotted line. 



This should not affect the comparison since the angular 
momentum profile does not change during Si-burning. 
One can see in our case that convergence only occurs 
above Vini=200kms^^. Indeed, the average specific an- 
gular momentum of the core (assuming a 1.7 Mq core) 
at the end of the calculation is 1.326 10^^, 1.801 10^^ and 
2.106 lO^^cm^s-i for Uini= 100, 200 and 300 kms^^ re- 
spectively. 

Concerning the evolution of the angular momentum, 
the general picture is the following. Mass loss removes an- 
gular momentum from the surface and transport processes 
(convection and rotational mixings) redistribute angular 
momentum inside the star (see Sect. 14.21 and the refer- 
ences given above for details). Here we are only concerned 
about the evolution of the angular momentum of the core 



of the star. During H-burning, both our models and mod- 
els without the inhibiting effect of the /^-gradient on mix- 
ing (models without "B") from HLWOO show a large de- 
crease of the angular momentum of the core. On the other 
hand, in HLWOO models including the inhibiting effect of 
the /^.-gradient (models with "B" ) , the core does not lose 
much angular momentum during H-burning. In our mod- 
els, thermal turbulence is taken into account and is able to 
overcome the inhibiting effect of the /i-gradient. HLWOO 
do not include the thermal effects and in their situation, 
the inhibiting effect of the /^-gradient is almost complete 
even with a reduction parameter /^ equal to 0.05. The dif- 
ferent treatment of rotation (and especially the different 
way the inhibiting effect of the /i-gradient is included) has 
therefore a strong impact on the evolution of the angular 
momentum of the core during the MS and explains the 
difference between the various models. 

At the end of H-burning, the core contracts and the en- 
velope expands. This restructuring phase is accompanied 
by a formation of a very deep external convective zone. 
At the same time, shell H-burning creates a short-lived 
intermediate convective zone. These changes may affect 
the angular momentum profile. The largest change in our 
models is the creation of a large drop of the angular mo- 
mentum at the bottom of the external convective zone (see 
Fig. [Tj) . This is due to the fact that convection enforces 
solid body rotation and therefore angular momentum is 
transported at the outer edge of the convective zone. No 
significant change is seen in the core. 

During He-burning, the trend is the opposite from H- 
burning. In both our models and models without the in- 
hibiting effect of the /t-gradient on mixing from HLWOO, 
the angular momentum in the core decreases slightly. On 
the other hand, in HLWOO models with the inhibiting ef- 
fect of the /i-gradient, the core loses a significant amount 
of angular momentum after H-burning. The reason is the 
following. During H-burning, in HLWOO models with the 
/i-gradient effects on mixing, even though the core does 
not lose much angular momentum, the layers just above 
it lose angular momentum (due to various transport pro- 
cesses). This creates a large angular velocity gradient at 
the edge of the core which increases rotational mixing 
during He-burning. Furthermore the successive convective 
and semiconvective zones (due to the restructuring phase 
and shell H-burning) mix as well the angular momentum 
of the outer parts of the core with layers above the core 
and a large amount of angular momentum is transfered out 
of the core at this time. In our models (as well as those 
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from HLWOO without /x-gradient), angular momentum is 
transfered to the layers above the core during H-burning. 
This creates a smaller gradient of angular velocity at the 
edge of the core at the end of H-burning and thus ro- 
tational mixing is weaker during He-burning. Therefore, 
the angular momentum of the core does not change as 
much during He-burning. Note also that He-burning is 
ten times shorter than H-burning and that there is less 
time to mix. As said in Sect. 14.21 during the advanced 
stages, the angular momentum profile does not change 
substantially. Only convective zones create spikes along 
the profile. 

The comparison of the final angular momentum profile 
of the different models, all with the same initial mass and 
surface angular velocity, is shown in Fig. [23 The thick 
solid line (HMM) corres ponds to our model. The models 
from lHeeer et all l|2000(l are drawn with a dotted-dashed 
line for model E20 (no /i-barrier) and with a dashed line 
for model E 20B (fx-barrier with /^ = 0.05). Finally, model 
m20b5 
magnetic 

the dotted line. Even though the evolution of angular mo- 
mentum differs between our model an d model E20B (wit h 
the /i-gradient effects on mixing) from Hcgc r et al.l (|2000() , 
the final value of the angular momentum of the core is very 
similar for these two models. This confirms the possibility 
of the for mation of GRBs v ia collap sars from rotating mas- 
sive stars l|Wooslev fc Heger.2003:.Heger et al.-2003.) if the 
effects of magnetic field (not included in our work) are 
small. Indeed, for example, the 25 Mq model ends up as 
a WR star with a core having enough angular momentum 
^ 10^^ c?7i^ s^^ for the material just outside the core, 
see Fig. O to create a collapsar. The difference between 
our model and r nodel E20 (without the /i-gradient effects 
on mixing) from lHeger et al. ( 2000) is probably due to the 
combination of the non-inclusion of the /i-gradient effects 
on mixing and of the different treatment of meridional cir- 
culation (see Sect. 17. 4. 2|) . From model m20b5 (see Fig.OBjI. 
one sees that the inc l usion of the effects of magnetic fields 
according to ISpruifi l)2002r) decreases significantly the fi- 
nal angular momentum of the core. In this situation, the 
core rotates too slowly and cannot produce a collapsar. 
We can also compare our models with the observed rota- 
tion period of young pulsars. Rotating models without the 
effects of magnetic fields have about 100 times more an- 
gular momentum at the pre-SN stage than the observed 

I "— -11 -1 

young pulsars IjHeger et al J 1200011 . M odels includin g the 

effects of magnetic fields according to Soruig 1)200211 have 
about 5-10 times more angular mome ntum at the pre-SN 
stage than the observed young pulsars l|Heger et al.l2003f ). 
This means that in any case, add itional slow dow n is nec- 
essary during th e core collapse IjWooslev fc Heger 200,Si 
iFrver fc Warrenl l2004') in order to reproduce the observed 
rotation periods of young pulsars. The question that needs 
to be answered is when and how this slow down occurs. 
Further developments will therefore be of great impor- 
tance for the formation of both NSs and GRBs. The topic 
of the final angular momentum of our models and its impli- 



cations for further evolution will be developed in a future 
article. 



7.6. Lower mass limit for models to reach iron core 
collapse 

As said above, our 12 Mq models have not been pur- 
sued beyond the O and Ne-burning phases for the rotat- 
ing and non-rotating models respectively. Nevertheless, 
we think that the rotating model has the potential to 
reach an iron core while the non-rotating rnodel does not. 
Recent calculations done bv lHeger et al.l l|2000ll are sim- 
ilar to ours on that point: Their non-rotating models as 
well as the rotating models E12B, F12B and G12B neither 
reach core-collapse. Only the model E12 reaches core col- 
lapse but the physics used in that last model does not 
include /i-gradient inhibiting effects on rotationally in- 
duced mixing. At the same time, rec ent n on-rotating 13 
Mq models in IWooslev et al.l ij2002l) and iLimongi et alJ 
pOOO) reach core-collapse. Therefore we expect the lower 
mass limit for non-rotating models to reach the standard 
iron core collapse to be around 12-13 Mq in agreement 
with Nomot o fc Hashimo to (1988). Our rotating models 
tend to show that this limit should be lower for rotating 
stars. A finer grid of models around 12 Mq would help 
constraining this limit. 



8. Conclusion 

The Geneva evolution code has been improved in order 
to model the pre-supernova evolution of rotating massive 
stars. We extended the nuclear reaction network with a 
multiple-a elements chain between carbon and nickel for 
the advanced burning stages. We also stabilized the in- 
ternal structure equations using Sugimoto's prescription 
( Sugimoto 1970 ) . Finally, we added dynamical shear to 
the other rotationally induced mixing processes (secular 
shear and meridional circulation). 

We calculated a grid of stellar models at solar metal- 
licity with and without rotation and with masses equal to 
12, 15, 20, 25, 40 and 60 Mq. 

Concerning the evolution of rotation itself during the 
advanced stages, the angular velocity increases regularly 
with the successive contraction of the core while the an- 
gular momentum does not change significantly (only con- 
vection creates spikes along its profile). This means that 
we can have a good estimate of the pre-coUapse angu- 
lar momentum at the end of He-burning. Gomparing our 
pre-SN models with the criteria for collapsar progenitors 
l)Wooslev fc HegeiJ l2003V we find that WR stars are pos- 
sible progenitors of collapsars. However, in this work we 
neglected the effects of magnetic fields. Further develop- 
ments will be very interesting for the formation of both 
GRBs and neutron stars. Dynamical shear, although very 
efficient, only smoothens sharp angular velocity gradients 
but does not transport angular momentum over great dis- 
tances. 
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We find that rotation significantly affects the pre- 
supernova models by the impact it has during H and He- 
burnings. We clearly see the two mass groups where either 
rotationally induced mixing dominates for M < 30Mq or 
rotationally increased mass loss dominates f or M > 30Mq 
as already discussed in iMevnet fc Maeden l)2003J) . 

We show that rotation affects the lower mass limits for 
the presence of convection during central carbon burning, 
for iron core collapse supernovae and for black hole forma- 
tion. The effects of rotation on pre-supernova models are 
most spectacular for stars between 15 and 25 Mq. Indeed, 
rotation changes the supernova type (lib or lb instead 
of II), the total size of progenitors (Blue instead of Red 
SuperGiant) and the core sizes by a factor ~ 1.5 (bigger in 
rotating models). For Wolf-Rayet stars (M > 30Mq) even 
if the pre-supernova models are not different between ro- 
tating and n on-rotating models, thei r previous evolution 
is different IjMevnet fc Maeden 1200,31) . We also compare 
our results with the literature. The biggest differences are 
the final mass and the various core masses. We obtain big- 
ger core masses and this should have a strong impact on 
yields. 

Future developments are planned to be able to follow 
the evolution until core collapse as well as follow neutron 
excess and detailed nucleosynthesis during the entire evo- 
lution. 
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